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Abstract. We investigate theoretically the possibility to observe dynamical mode 
locking, in the form of Shapiro steps, when a time-periodic potential or force 
modulation is applied to a two-dimensional (2D) lattice of colloidal particles that are 
dragged by an external force over an optically generated periodic potential. Here 
we present realistic molecular dynamics simulations of a 2D experimental setup, 
where the colloid sliding is realized through the motion of soliton lines between 
locally commensurate patches or domains, and where the Shapiro steps are predicted 
and analyzed. Interestingly, the jump between one step and the next is seen to 
correspond to a fixed number of colloids jumping from one patch to the next, across 
the soliton line boundary, during each AC cycle. In addition to ordinary “integer” 
steps, coinciding here with the synchronous rigid advancement of the whole colloid 
monolayer, our main prediction is the existence of additional smaller “subharmonic” 
steps due to localized solitonic regions of incommensurate layers executing synchronized 
slips, while the majority of the colloids remains pinned to a potential minimum. 
The current availability and wide parameter tunability of colloid monolayers makes 
these predictions potentially easy to access in an experimentally rich 2D geometrical 
configuration. 


PACS numbers: 82.70.Dd,83.85.Vb,68.35.Af 
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1. Introduction 

Driven nonlinear systems can show widely intrigning patterns of varied and often 
nnexpected dynamical behavior. An especially interesting case is that of many identical 
interacting and crystallized particles moving in a periodic external potential: a common 
sitnation, e.g. for adatom strnctnres sliding on a crystalline snrface. In the special 
case when, besides the main time-independent force which canses the sliding, an 
additional external ”AC” force component is present and oscillates periodically in time, 
synchronization phenomena are known to occnr, depending on the relative importance 
of the interactions and on the geometrical arrangement of the particles. 

Synchronization phenomena of this kind have long been described nsing simple, low¬ 
dimensional, phenomenological models, which often snffice to captnre the main featnres 
of the complex dynamics involved. The basic example is provided by the simple yet 
fnndamental one-dimensional (ID) Frenkel-Kontorova (FK) chain model [Il|2l[3] whose 
soliton, or kink pattern is known to exhibit a nontrivial intermittent dynamics, when 
additional time-periodic forces act on the otherwise steadily sliding chain |1] . This simple 
ID model however is not generally snfficient to describe in fnll the complex dynamics 
of atoms at the interface of two materials in relative motion. Experimental stndies 
of two-dimensional (2D) arrays of mntnally repelling colloids pnshed across a periodic 
“corrngation” potential landscape have recently been carried ont, providing a closer 
analog of the sliding of crystalline interfaces [SlEllTllH]. Before sliding, the observed 2D 
patterns of misfit dislocation lines (referred to as solitons or kinks) are very mnch akin 
to the Moire patterns formed by atomic overlayers over crystal snrfaces. Besides their 
snrface physics analog, the colloid systems have an interest on their own becanse they 
provide a ready parameter-controlled system for the stndy of friction between crystalline 
snrfaces iiiTiiig. 

Here we explore, ahead of experiments, jnst very recently available in a simpler ID 
geometry [H], a sitnation where synchronization, whence the existence and the actnal 
dynamical natnre of Shapiro steps, is songht by oscillating in time either the pnshing 
force or the amplitnde of the {x, y) periodic corrngation potential. In particniar we focns 
on the second, corresponding to an ideal setnp where colloids thread a partially time- 
modnlated spatially-periodic optical potential, and investigate in particniar the effects 
of the mismatch between the colloid-colloid average spacing of the 2D colloid crystal 
and the periodicity of the optical potential. 

Shapiro steps were initially observed [12] in the context of Josephson jnnctions. 
The crncial ingredient at the start was Josephson’s prediction [13] of a coherent 
cnrrent Asin^ flowing between two snpercondnctors separated by a thin insnlating 
barrier in presence of a phase difference (f) = 4>i — 4>2 between the two snpercondnctor 
order parameters. Moreover, a potential difference ^ was predicted to be 

associated to a phase 0 that changes in time, with the remarkable conseqnence, known 
as AC Josephson effect, that the coherent cnrrent wonld oscillate with a freqnency 
Uq = 2eVdc/h in presence of a constant DC bias Vdc- The fingerprint of ojq can be 
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detected through an ingenious resonance experiment [12], in which a current I{t) = 
hdc + -^acsin(a;t) is passed through the junction. A simple model jlH [I5| capturing 
resistive and capacitive effects in the junction would allow us to write: 


hC df(j) h d(j) 
2e df^ 2eR dt 


+ Ic sin f = I(t), 


( 1 ) 


where C is the capacitance of the junction, and the corresponding current, 

while ^ accounts for the resistive part of the current, due to quasi-particles. 

In general, a fit) = uj^t + A0(t) where 2S.^{t) is a bounded function, implies a (time- 
averaged) DC voltage Vdc = ^ across the junction. In a non-resonant situation, sin fit) 
oscillates in an aperiodic way with zero average, (sin0(t))time = 0; hence the (time- 
averaged) DC voltage predicted by Eq. ([1]) follows the usual Ohmic I-V characteristics, 
Vdc = Ride- Resonances instead can make (sin0(f))time 7^ 0 by the following phase¬ 
locking mechanism US]: the bounded variation A0(t) becomes a periodic function of 
time, characterized by the frequency u of the external driving, while the slope of 
0(f) — the frequency ujq — exactly matches a multiple of the driving frequency cu: 
(Uq = nuj. In such resonant situations, sin 0(f) becomes periodic and with non-zero 
average, (sin 0(f))time 7^ 0, implying that the constant terms in Eq. ([T]) obtained by 
averaging over a period sum up to: 


„ “1“ A(sin 0(f))tijcie hdc ■ (2) 

H 

In essence, at resonances, Jdc can be changed by an amount ~ A (sin 0(f)) time without 
affecting the voltage bias Vdc: these are the Shapiro steps in the I-V characteristics, 
where the the voltage step is strictly quantized through a finite range of current. 

A mechanical analogy of the Josephson junction model of Eq. ([1]) is straightforward, 
with the phase difference 0(f) playing the role of a coordinate x{t), the voltage V oc d(()/dt 
that of the velocity of a particle moving in a washboard potential, and the current 
/(f) the role of an external force. Shapiro steps have then a natural translation into 
“quantized” steps of the velocity of the moving particle for a finite range of applied force, 
when the washboard frequency goes in resonance with the AC perturbing frequency u. 
This resonance phenomenon is robust, and one could for instance periodically modulate 
the amplitude of the space-periodic corrugation potential [the analog of R in Eq. ([1])] 
rather than the external force, and the same physics would follow. 

As is shown in Eq. (IT|), the Shapiro physics involves a single degree of 
freedom. However, in systems of many interacting particles, crystallized but not rigid, 
synchronization phenomena of the Shapiro kind are likely to give rise to additional 
nontrivial dynamical effects due to concerted multi-particle motion. The strong interest 
in synchronization phenomena in multiple physical contexts has generated numerous 
reseach works in recent years, among which we recall Refs. mi HHEHllIHl EHl I2I11221 
l23l l2n |25] . The system where we propose to explore the Shapiro physics consists of a 
collection of repulsive colloidal particles dragged by an external force over a spatially 
periodic corrugation potential, whose amplitude includes both a static and a time- 
oscillating part. When the frequency u of the oscillation is a multiple of the characteristic 
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(washboard) frequency Uq, the forced sliding motion of the 2D colloid lattice can give rise 
to Shapiro steps. Combined experimental and simulation studies already addressed in 
the past the effect of synchronization and the ensuing Shapiro steps in Brownian particle 
dynamics [23]. Other simulation work [26] addressed the case where the 2D sliding 
lattice and the periodic potential are fully matched, in addition to ratcheting conditions 
in mixtures. Very recently, the microscopic dynamics underlying mode locking in a 
colloidal model system has been recorded in a simple, yet noteworthy, ID experimental 
setup [H]. In this work we study the feasibility of observing Shapiro steps under 
reasonably realistic conditions in 2D sliding colloid monolayers, under the wider range 
of conditions that can be realized experimentally. More specihcally, we will describe a 
planar system of monodisperse repulsive particles, forming a 2D triangular lattice with 
spacing acoii, sliding within a periodic corrugation potential due to an optical lattice, also 
triangular with spacing oias, whose amplitude is partly modulated periodically in time, 
under various conditions as determined by different choices of the commensurability 
ratio p = aias/ocoii- Using classical molecular dynamic we will monitor the individual 
motion of the particles as well as that of their center of mass (CM) under the action of 
a DC force. As a function of this force the forward CM velocity of the colloidal particles 
should, owing to Shapiro synchronization, become ’quantized’ in step-like structures as 
a function of the applied force. Our aim is therefore to seek, detect and characterize 
these step structures in relation to the absence (in the fully matched, commensurate 
case p = 1) or presence, e.g., for p < 1 of pre-existing soliton-like defects or kinks due 
to mismatched relative spacings of the colloidal lattice and of the periodic potential. 

2. The model 

Molecular dynamics (MD) simulations are based on the same model already introduced 
[To] to describe the motion of mutually repulsive charged colloidal particles driven across 
a periodic potential generated by a light interference pattern, as realized in experiments 
by Bohlein et al. [ZllH]. In short, we describe the charged colloidal particles as classical 
point-like objects, whose dynamics is affected by their mutual repulsion, the action of 
external forces, plus the interaction with the viscous fluid where they are immersed. 

The equation of motion for the j-th particle is 

mrj{t) + virjit) - Vdx) = -Vrj(U2 + Uext) + fj(t ), (3) 

where rj is the position of the j-th colloid in the 2D plane where it can move, rj is 
the friction coefficient determined by the effective viscosity of the fluid in which the 
colloids are immersed, and Vd is the fluid drift velocity, giving rise to a Stokes driving 
force F = rjVd, pushing forward all colloidal particles. U 2 is the two-body inter-particle 
potential; Uext is the corrugation potential describing the interaction with the spatially 
periodic light-held pattern. In this work we especially focus on the effect of time- 
periodic oscillations of the amplitude of Uext- The hnite-temperature Brownian motion 
of colloids is introduced in a standard Langevin approach, involving the viscous friction 
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Table 1: Adopted simulation parameters, expressed in model units |T0]. For example, 
the simulation-cell sides and Ly are given in units of the average colloid lattice spacing 
flcoii- A^minima IS the number of identical local minima of the corrugation function W{y) 
contained in a simulation cell. 


term r]{rj{t) — u^x), plus the Gaussian random force ij(t) [27j. Due to the low overall 
velocity {vd — 10 /ims“^) of the colloidal particles, the inertial term can be neglected, 
and an overdamped diffusive motion is reasonably assumed, with a sufficiently large rj. 
In our work rj = 28 expressed in the same model units described in Ref. [10]. Typical 
simulation parameters are reported in Table [H 
The two-body interaction potential is 

N 

EU„k(|rj-r,,|), (4) 

j<j' 

where the screened Coulomb repulsion R is a Yukawa potential: 

Wuk(’^) = —exp(-r/AD) • (5) 

r 

The average nearest-neighbor separation of colloids in Ref. [7] is r = Ucoii — 5.7 pm ~ 
SOAi). Under conhnement, this repulsion establishes a triangular lattice of colloids, 
which can be thought as reasonably defect free at temperatures not too high and sizes 
sufficiently small to render the possibility of Nelson-Halperin dislocations [2S] irrelevant 
to the present case. 

The one-body term 

N 

Uext = , (6) 

j 

describes the interaction of the colloids with the 2D spatially periodic potential, 
representing the optical lattice corrugation. Its actual form, determined by Rxt(r), 
could in principle be shaped with some freedom, but the simplest sinusoidal form is 
sufficient to describe the main phenomena. We thus take: 


Uext(r, t) = [Vo + Ao sin((:ut)] W (r), 


(7) 


where 


»'(r) = -- 


3-1-4 cos 


27rr,, 


cos 


27rra 

®las 
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47rr,, 




®las / 


( 8 ) 


. \/3 Ulas 

a unit-amplitude eggcarton-type periodic potential of 6-fold symmetry representing the 
optical-lattice triangular corrugation. The time-dependent amplitude Vq -|- Aosin(c(;t) 
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gives to the corrugation an additional sinusoidal time dependence that can induce 
synchronization effects. Its amplitude 2Ao yields a total corrugation ranging between 
Vo — Aq and Vq + Aq with an AC modulation frequency v = u/{2tt). The minimum 
force required to dislodge an isolated colloid at T = 0 along direction x is at time t 
Svr Vo + Aq sin(a;f) 

9 ®las 


Fsi = 


(9) 


The quantity 8Vb/(9aias) is used as our unit of force. In general, the spacing aias of the 
corrugation lattice potential W{x, y) can be tuned in such a way to be either matched 
or mismatched with the colloidal lattice. In mismatched configurations, one observes, 
both in simulations [10] and in experiments [7] , the formation of misfit dislocations also 
called topological solitons, or kinks/antikinks in the language of the FK model [29] • Their 
existence and their motion under the external force dominates the frictional properties 
of the sliding lattice, since it provides the mechanism for mass transport [30] . 

As for the oscillation frequency, it needs to be low enough for the overdamped 
colloidal system to be allowed enough time to follow adiabatically the AC modulation 
in the corrugation-amplitude oscillation. The specific condition is z/ 77/m, and we use 
z/ < 0.1 in model units. Shapiro-step structures can arise out of synchronization of the 
washboard frequency ujq of the colloids driven by F over the corrugation potential, with 
the time-modulation frequency u of Eq. (0. At finite temperature (e.g. in experimental 
conditions), v must also not to be so small that particle random diffusion overshadows 
the ordered synchronized crossing of potential energy barrier. 

Our protocol begins with a preliminary relaxation preparing the sample in its 
equilibrium, force-free initial state. Sliding simulations are then conducted by applying 
a force F = pvd along x on each particle, similar to the viscous drag of the flowing fluid 
in experiments. The static force F is slowly ramped upwards in small steps AF, so as 
to explore the desired force range. At each force value, after an initial transient lasting 
between 1 and 10 oscillation periods, the particle CM displacement Axcm is extracted 
for the successive simulation time r, typically a large integer (of order 100) number 
of periods z/“^ of the time-dependent perturbation. We then record the average speed 
Vera = Axem/'T for each value of the driving force F. 

Experimental conditions [7] actually involve an inhomogeneous 2D configuration, 
whereby only a circular central portion of the colloid monolayer raft is submitted to 
the laser field with its associated optical lattice potential. However, only phenomena in 
the central part of this circle are eventually studied, and there the colloids are virtually 
homogeneous. Therefore, we conduct the present study in a fully homogeneous 2D 
colloid lattice with periodic boundary conditions (PBC) qualitatively represented in 
Fig. [H Given a cell size and a particle number, we will simulate either the commensurate 
and fully matched geometry (p = 1), or mismatched underdense geometries (p < 1), by 
tuning appropriately the corrugation lattice spacing oias relative to that of the colloid 
lattice, Ocoii- 
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Figure 1: The x Ly = ISocoii x 15\/3acoii rectangular simulation cell (solid line) 
adopted for the p = 3/4 simulations. In this initial snapshot, the 450 particles still form 
a perfect triangular lattice. 


3. Perfect lattice matching 

We begin with the fully matched commensurate geometry where, at p = 1, the two 
triangular lattices coincide, and no preliminary relaxation is needed, each particle falling 
in a potential minimum. We then turn on the sliding force with a time-dependent 
modulation of amplitude Aq = 40% Vo, with an intermediate corrugation strength 
Vo = 0.5, at T = 0. The average CM velocity as a function of the external force 
F is presented in Fig. O As a function of the sliding force, it shows Shapiro steps, 
i.e. plateaus in the average CM velocity. These steps correspond to integer multiples 
of z/oias in the force-velocity response. In these T = 0 simulations the motion of the 
colloids is deterministic, leading to perfectly flat steps with null error bars, as long as 
the CM speed is averaged over an integer number of periods The reference velocity 
naias is that of a particle advancing by one corrugation spacing aias in one oscillation 
period Depending on the applied DC force, all particles advance together by an 
integer number of lattice spacings aias at each period of the AC modulation. As a 
function of n both the velocity jump between successive Shapiro steps, and the force 
step width increase proportionally to u, as expected. Correspondingly, the number of 
steps in a given force range decreases for increasing z/. In the large force regime, for 
^ (Yo + Ao)/aias — 1.95, the steps rapidly decrease in width, and merge into a 
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Figure 2: T = 0 Shapiro steps observed in the CM velocity parallel to the applied force 
obtained for a fully commensurate colloid lattice (aias = acoii, P = 1) with a corrugation 
amplitude Vq = 0.5, a modulation amplitude Aq = 0.2, for values of the frequency 
n = 0.01, 0.02, and 0.03 (in model units). The velocity is scaled by z/aias, to emphasize 
the proportionality of the exactly quantized steps to the modulation frequency. 


smooth velocity increase proportional to force, characteristic of a viscous regime, similar 
to the non-modulated Aq = 0 dynamics mm- 

In this zero-temperature lattice-matched system, all colloids execute the same 
movement to negotiate the crossing of the same energy barriers. As a result, the 
colloid-colloid mutual spacing remains constant, and the colloid-colloid internal force 
does no work. The dynamics would be exactly the same if the colloidal particles were 
not interacting, or if the simulation involved one particle only, moving in the same 
potential-energy prohle. As a check, we repeated the simulations with a single particle, 
indeed recovering exactly the same patterns as in Fig. [2l 

We explore now how the AC modulation amplitude Aq affects these Shapiro steps. 
Figure [3] shows the velocity curves obtained with the same average corrugation amplitude 
(Vq = 0.5) and AC modulation amplitude rising from the same initial value Aq = 40% Vq 
as in the central panel of Fig. |2] (z/ = 0.02), to larger values Aq = 70% Vq and 
Aq = 100% Vo, at the same frequency. For larger Aq the colloid system starts to move 
earlier, i.e., at a smaller DC force than for a smaller Aq. Note that at 100% relative 
modulation amplitude, every oscillation period has a brief instant where the colloid 
layer actually slides freely. Despite this free-sliding instant, the force needed for the 
monolayer to depin and reach the hrst Shapiro step, even though smaller than at 40% 
modulation amplitude, appears anyway to be nonzero. The reason for that is the hnite 
time it takes the colloid to drift across the momentarily turned off barrier to the next 
substrate minimum, when the driving force F is small. This issue is discussed in detail 
in Sect. |5] below. 
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Figure 3: Shapiro steps for z/ = 0.02, Vq = 0.5, for increasing AC modulation amplitude 
Aq. Pluses: Aq = 40% Vq, same as in the central panel of Fig. [2l Crosses: Aq = 70% Vq. 
Squares: Aq = 100% Vq. The step heights are unchanged, their widths increase, and the 
jump between successive steps becomes steeper for increasing Aq. 



Figure 4: Finite-temperature Shapiro-step structure for a fully commensurate system. 
Temperature T = T^oom, and the other simulation conditions {u = 0.02, Vq = 0.5, 
Aq = 70% Vo) are the same as for the crosses data of Fig. |3l The effect of thermal 
motion is barely detectable, apart from a slight rounding on the step edges. Thermal 
fluctuations of the CM velocity, reflected by the error bars, are visible at rises between 
steps, but negligible inside the steps. 


3.1. Thermal effects 

So far we worked at T = 0, because that regime yields a clearer picture. Next, however, 
we can ask how the Shapiro steps might be affected by thermal fluctuations in the 
experimental, room-temperature conditions. To address this question, we perform new 
simulations at room temperature T = T^oom (Troom = 0.04 in system units), still at 
p = 1 and keeping all other parameters (mimicking the experimental colloid system of 
Ref. [7]) the same. Figure H] shows that the Shapiro-step structure is barely affected; 
Troom is a very moderate temperature for this system. Thermal fluctuations (estimated 
by the standard deviation of Ucm obtained over several sliding simulations realized with 
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independent random-number sequences and different initial configurations ensuring a 
better sampling of the thermal equilibrium state) are visible, as reflected by the error 
bars, at jumps between steps. Thermal fluctuations produce a rounding of the plateau 
onset, but are barely visible inside each plateau. 

4. Mismatched lattice spacing 

When the commensurability ratio p = aias/flcoii differs from unity the colloid and 
the corrugation lattices are mismatched. As we shall see, in this case the internal 
structure of the sliding lattice plays a role, and the expected Shapiro-step phenomena 
become more interesting. It must be clarihed at the outset that the Shapiro steps 
exist only if the system is pinned by a finite static friction force. For example, a 
2D colloid lattice incommensurate with a weak periodic corrugation is an unpinned 
system which can slide “superlubrically” without static friction, and will exhibit no 
Shapiro steps. The same incommensurate lattice must however become pinned when 
the corrugation amplitude exceeds a hrst-order depinning-pinning transition threshold 
value [31], and here Shapiro steps can exist. This is the situation which we concentrate 
upon here. We describe in detail two different examples of mismatched configurations, 
p = 14/15 ~ 0.93 and p = 3/4 = 0.75. Both are underdense systems compared to the 
perfectly matched one with p = 1 considered in Sec. [S] and both are of course rationally 
commensurate. However, the commensurabilty of p = 3/4 is so to speak stronger, that 
in p = 14/15 weaker, the latter case in a sense closer to incommensurability, itself a 
condition impossible to reach in a finite-size PBC realization. We will further assume 
that the colloids and the corrugation lattice do not undergo a relative rotation, e.g., 
of the type studied in Refs. [32l [33l [Ml iMl [36], and therefore remain aligned. With 
that stipulation we adopt for p = 14/15 a rectangular periodic simulation cell involving 
392 colloids threading 450 corrugation potential minima. With 14 colloidal particles 
threading 15 potential wells along each principal direction, we obtain a widely spaced 
hexagonal network of misfit dislocation lines {antisolitons), where vacancies segregate. 
Details of this case will be shown later. For p = 3/4 we also adopt a rectangular 
supercell, here involving 450 colloids and 800 potential minima, shown in Fig. [T] The 
higher vacancy concentration leads to a greater mismatch and therefore a denser packing 
of antisolitons. 

In these mismatched systems we must first of all identify the corrugation amplitude 
at which, as was mentioned above, a pinning/depinning transition takes place in the 
absence of AC modulation (Aq = 0). In the limit of inhnitesimal applied force, 
this transition was recently discovered and characterized by Mandelli et al. |31j . 
Ignoring here the weak 14/15 commensurability and taking this to represent a truly 
incommensurate case, this is the 2D analog of the celebrated ID “Aubry transition” 
[29] EHl El] here taking place at Vocrit — 0.4 [39]. With this stipulation, for every 
corrugation value Vq < focrit the colloid sliding over the incommensurate corrugation is 
“superlubric”, and will take place for arbitrarily small applied force. By contrast, when 
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Figure 5: DC Force F - corrugation amplitude Vq phase diagram for p = 14/15 
mismatched colloids, in the absence of AC modulation (Aq = 0). The red dots and line 
represent the depinning line where at T = 0 the colloid abandons the statically pinned 
state (shaded yellow region) in favor of a (nearly) free-sliding state for larger force F. 
A truly incommensurate monolayer, where p is irrational, would exhibit a very similar 
transition line, the 2D analogue of the ID Aubry transition [311 [37] . Horizontal solid and 
dot-dashed lines: a choice of corrugation (Vq = 0.5) and AC modulation Aq = 15% Vq 
amplitudes. In the force range for which both dot-dashed lines are in the shaded region, 
the colloid layer is pinned. For large force, when both dot-dashed lines lie below the 
depinning curve, the colloids advance freely. The intermediate force interval is the one 
affected by the synchronization to the AC potential modulation, giving rise to Shapiro 
steps. 


Vo > Vbcrit, the monolayer is pinned, and a hnite static friction force F > 0 is required 
to make the system slide [3T]. The curve of Fig. |5] represents precisely this depinning 
transition line for the p = 14/15 mismatched system, for nonzero and increasing force. 

To understand what happens when we superpose an AC corrugation modulation 
Aq > 0 to the “phase diagram” of Fig.jSl we add horizontal dot-dashed lines representing 
the maximum and minimum value covered by the overall corrugation magnitude during 
an oscillation. As an example, in Fig. [5] we sketch Vq = 0.5 (horizontal solid line) and 
Aq = 15% Vo (horizontal dot-dashed lines). Synchronization between the washboard 
sliding frequency and the AC modulation, and consequently Shapiro steps, can arise 
naturally in the intermediate range of forces 0.02 < F < 0.15, where at each modulation 
period the system crosses twice the pinned-unpinned transition curve. Figure [5] shows 
that depending on Vq and Aq, the F dependence can vary quite substantially. For 
example, if Vq -|- Aq < Vocut, then for any F the colloids will slide freely. For 
%o +Ao > Vocrit but Vq — Aq < Vocrit the initial pinned region (F < 0.02 in hgure) 
would disappear: in this case the first Shapiro step could extend to arbitrarily small 
applied force, by reducing the modulation frequency z/, as further discussed in Sec. |5l 
This phase diagram can be taken as a guide for the choice of the simulation parameters. 
In order to observe Shapiro steps, the colloid lattice must be pinned, e.g., Vq must be 
above the critical corrugation Vocrit- The range of force where Shapiro steps exist widens 
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Figure 6: The Shapiro-step structure of the lattice mismatched p = 14/15 conhguration 
simulated at T = 0, with corrugation amplitude Vq = 1, modulation Aq = 60% Vq, and 
AC frequency z/ = 0.02. (a) A zoom-in of the region of the “trivial” commensurate 
step Ucm/(^aias) = 1 Corresponding to all colloids advancing by aias at each period, 
(b) A zoom-in of the region of the subharmonic steps, highlighting those matching 
w = |1 — p~^\ = 1/14 and integer submultiples thereof. 


out as the modulation amplitude Aq is increased. The phase diagram for p = 3/4 is 
quite similar to that of p = 14/15, and in particular Vo,crit — 0.42 is nearly the same. 

We now explore and describe the existence and nature of Shapiro steps in these 
mismatched cases, adopting an arbitrary but reasonable corrugation Vq = 1 > Vo,crit 
and AC modulation Aq = 60% Vq (different parameters were explored too). Unlike the 
p = 1 case, here a preliminary and careful relaxation of the initial force-free, equilibrium 
structure is required. An initial simulation for a duration r = 1000 is thus performed 
with a time-independent corrugation (Vq turned on, Aq = 0), and no external force 
applied {F = 0). Starting from the relaxed conhguration, a sequence of force-driven 
sliding simulations is run, now with the AC corrugation modulation turned on. 

4.1. p= 14/15 

4 . 1 . 1 . Step structure. Let ns consider p < 1, where it is useful to introduce as a misht 
measure the quantity 



®coll 


, 1 

w = 

1- 




O'las 


/ 


This quantity was found earlier to be the relevant parameter in the context of simulations 
of mismatched crystalline systems in mutual contact. Refs. |10l SH HJl HHl HU HSl 1461 
HU SHI sni EDI El]. The results of sliding simulation for w = 1/14, corresponding to 
p = 14/15, are shown in Fig. El The Shapiro steps exist here too, but the mismatched 
geometry exhibits a far richer pattern of Shapiro steps than the lattice-matched case: 
we detect tens of steps in a much narrower force range than in the commensurate case 
of Fig. m 
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First off, the regular integer plateau at Ucm/(^ ^las) = 1, at the same speed as the 
hrst step in fully matched geometry, Fig. |3l is detailed in the inset Fig. |6^. Additionally, 
Fig. [6] exhibits several other subharmonic steps at smaller driving force. It has long been 
known that for mismatched lattices also fractional subharmonic Shapiro plateaus, i.e. 
plateau velocities which are not integer multiples of v aias [H [52] are to be expected. In 
particular the inset Fig. [6 ]d zooms around a much lower speed, given by 


^ ®las 

Here tc = 1/14 and velocity steps appear at wm/n for several m and n values. All these 
Vera < ^flias plateaus are subharmonic Shapiro steps. Unlike in the integer plateaus, 
where all N particles advance by one (or several) lattice spacing(s) oias in one period 
at a fractional plateau characterized by Vera/ (j^nias) = wm/n, a total fraction mwN 
of all colloids advances by one lattice spacings during n modulation periods. 

The quantization of sliding velocity to the value of Eq. ca was discovered and 
studied in the context of the sliding of mismatched crystalline systems in mutual contact, 
and reported in the study of two sliding surfaces with an atomically thin solid lubricant 
layer in between [IHl HU SSI 1131 HU SSI SHI SIl SSI SSI IHDl ISI]- There, the ratio p 
was that of the lattice spacing of a “bottom” slider to that of the lubricant layer, w 
was the mean velocity at which the lubricant moved over the bottom slider, in units 
of the externally imposed speed of the “top” slider. The sliding of the lubricant with 
relative speed w arose in that case because the vacancies (or interstitials) forming the 
antisolitons (solitons) were the real entities dragged forward at the full speed of the top 
layer, all (other) particles remaining essentially static. In the following we shall clarify 
the relationship between the driven soliton dynamics in the sliding-friction models and 
the subharmonic steps in the dynamics of colloids driven in a modulated-amplitude 
potential. 

4.1.2. Subharmonically sliding antisolitons. Figure [7] shows that the antisoliton lines 
divide the cell into domains inside which the particles are essentially static at a minimum 
of the corrugation potential. The overall forward sliding motion demanded by the 
force F is actuated, for p < 1, through the backwards motion of the antisolitons, 
where all vacancies are lumped. The antisolitons move to the left in synchrony 
with the modulation, advancing faster when the instantaneous corrugation amplitude 
Vq + ^o sm{u}t) of Eq. ([7]) is minimum. At that moment a fraction of colloids jumps from 
the left to the right side of an antisoliton line, while all others remain essentially static. 
In the w m/n step, this colloid displacement results in the full antisoliton pattern moving 
by moias to the left every n modulation periods. Figure [8] illustrates the dynamics of a 
single colloid. Most of the time the particle just oscillates, following the AC modulation: 
when the amplitude Vq-I-Ao sm{ujt) is maximum, the particle is driven toward the nearest 
corrugation-potential well; when the amplitude is minimum, the particle relaxes due to 
the two-body repulsion. The amplitude of this oscillation is minimal when the colloid is 
near the center of a domain. Fig. [Hb, and increases as the colloid reaches the antisoliton. 
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Figure 7: p = 14/15, T = 0, tf/3 fractional plateau. The antisoliton pattern divides the 
colloid lattice into domains or patches inside which colloids are approximately matched 
to the static corrugation lattice. The underdense colloid monolayer slides rightward 
due to the leftward motion of the antisoliton line, while particles inside the domains do 
not slide. The specificity of this plateau is the behavior of particles at the antisoliton 
boundary. Left and right are two successive snapshots separated by a single oscillation 
period . At the end of this time interval, most colloids still occupy their initial 
positions, but the 12 highlighted particles (red) crossed the antisoliton line to join the 
next commensurate domain to the right. In the two following periods the remaining 
Nw — 12 = 16 particles will also cross, so that after the time nv~^ (here with n = 3) all 
Nw = 28 particles will have crossed. Different plateaus with different w and different 
subharmonic n will have an analogous but different pattern of traversing particles. 


which is the domain boundary. A positive, uncompensated velocity spike, see Fig. [HtL,b, 
signals the moment when the particle crosses the gap and enters the neighboring domain. 
Fig. It should be noted here that some fraction of the colloid rightward motion 

also takes place while inside the commensurate domain, which is not rigid and thus 
becomes progressively deformed elastically. 

We are now in a position to clarify the nature of the w subharmonic steps. To 
do this, we observe the colloidal pattern in a ’stroboscopic’ approach, by comparing 
snapshots taken at successive times separated by one period, namely when the 
corrugation amplitude, Eq. o. acquires its maximum value Vq + Aq, i.e. for successive 
half-odd integer values of ut. We highlight the particles that during the last period 
have underwent a displacement exceeding a certain threshold 6. Due to the elastic 
domain deformations discussed above this threshold must be chosen appropriately: for 
the conditions considered here a value 0.12acoii < ^ < 0.61acou provides consistent 
results for all plateaus. The highlighted particles represent colloids, whose number we 
designate as rui, which cross the antisoliton boundary between two neighboring domains 
in that period i of oscillation. In the fundamental plateau where the quantized speed 
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Figure 8: p = 14/15, T = 0: The time evolution of the x component of (a) 
velocity and (b) position of a single (highlighted in panels c-e) colloid, when F = 0.24 
and the CM velocity is in the fractional Shapiro step of Fig. |6l characterized by 
'?^cm/(i^aias) = w = 1/14. Because p < 1, the forward sliding motion is realized by a 
backward motion of vacancies, here forming antisoliton lines. Velocity shows alternating 
positive and negative spikes, plus a large positive peak when the leftward moving domain 
boundary reaches, between t = 465 and t = 515, the colloid, which then jumps across 
the antisoliton line, joining the next domain. Near the center of the domain - e.g. at 
t ~ 175, panel (c) - the lateral oscillation amplitude reaches its minimum, because the 
particle sits close to a potential well bottom. 


equals w, the pattern of advancing particles is the same at all oscillation periods, as 
seen, e.g., in Fig. [9] and equals Nw where N is the total number of particles. In the 
present case of Fig.|9]V = 392, w = 1/14, the Shapiro plateau has dimensionless velocity 
w, and a total of Nw = 28 particles cross the boundary to the next domain at each 
oscillation period. More specihcally, 22 colloids at the edge of the island aligned along 
y move straight in the x direction into the next corrugation minimum at a distance aias- 
In addition, 12 more particles along the oblique antisoliton line advance with a sidewise 
y component, thus proceeding into a minimum which is aias/2 to the right. The total 
movement amounts therefore to effectively 22 + 12/2 = 28 = Nw particles advancing 
by flias in each period. The very same picture applies to sliding a monolayer with any 
different lattice-spacing ratio p = (1 + The fundamental subharmonic step occurs 

again at velocity w. 

The pattern of advancing colloids is somewhat more complex in the other 
subharmonic steps. We explore these differences in Figs. [9] and [TOl To exemplify, let 
us consider again w = 1/14, but now focus on the w/2 plateau. Fig. [101 Here, a total 
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Figure 9: p = 14/15, N = 392 particles, Vq = 1, forward force F = 0.24, step 
tc = 1/14. Successive snapshot patterns of advancing colloids, separated by intervals of 
one oscillation period The antisolitons drift to the left, causing the CM to move to 
the right. The particles whose displacements exceed a hxed threshold 6 are highlighted 
(lighter and bigger). In this step 22 particles shift to the right (i.e., along x) in each 
oscillation period, crossing the vertical section of the antisoliton line, plus another 12 
cross an oblique section, with a lateral y component of motion, and only half distance 
along X. As a result a total of 22 + 12/2 = 28 particles advance by oias, joining the next 
commensurate domain. 


Time 



Figure 10: The subharmonic plateau w/2 for p = 14/15, w = 1/14, F = 0.185. 
Successive patterns of advancing colloids, separated by intervals of one oscillation period 
illustrating the propagation of an antisoliton in this Shapiro plateau. In a hrst 
displacement 10 + 4/2 = 12 (highlighted) particles jump (the 4 particles crossing oblique 
antisolitons advance by a half lattice spacing along x). Then, at the next period, another 
12 + 8/2 = 16 (highlighted) particles complete the advancement of the string of particles 
bordering the antisoliton line. At successive periods, the alternating advancement of 12 
and 16 particles repeats itself involving successive lines of colloids more and more to the 
left. 
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particle number Nw again jumps between a domain and the next: but that takes place 
in two successive oscillation periods instead of one as in the harmonic w step. In the 
hrst AC oscillation period 10 + 4/2 = 12 colloids cross the antisoliton line, whereas the 
remaining 12 + 8/2 = 16 jump in the second period. Similarly, in the step leading to the 
w/3 plateau, it takes a sequence of three successive oscillation periods to complete the 
advancement which transfers Nw particles across the antisoliton. The partial transfers 
are in this case 10 + 4/2 = 12, 8 + 4/2 = 10, and 4 + 4/2 = 6 particles in period 1, 2 
and 3 respectively. 

In summary, each Shapiro step is characterized by a periodic advancement pattern 
during which an integer multiple of Nw particles jump from a commensurate domain 
to the next. However the m/n-subharmonic step of quantized velocity 

Ucm rn 

- = w — 

^ ^las ^ 

is accomplished as a composite of n successive modulation periods, each lasting of 
the AC modulation. After one such complete migration, the entire pattern of positions 
is displaced bodily by m lattice spacings aias- The subharmonic step leading to wm/n 
is thus composite, involving n AC modulation periods. The total number mNw of 
particles that cross the antisoliton is partitioned between the n periods, 

n 

mNw = y/ . (13) 

i=\ 

The detailed numbers m* of crossing particles at each individual AC period i will 
generally depend on p, Vq, Aq, and on the angle between the sliding direction and 
the antisoliton (soliton) orientation. Occasionally, the partition flT^ may even change 
along a subharmonic step as F increases. We have not found a general analytic rule, if 
any, determining this partition. 



4.1.3. Effect of the corrugation and modulation amplitudes. We now investigate the 
role played by the corrugation amplitude Vq and by the modulation amplitude Aq in 
the formation of Shapiro-step structures as exemplihed in the p = 14/15 mismatched 
conhguration. First of all we remind ourselves that below the pinning threshold 
corrugation amplitude Vocnt (which, as seen on Fig. [5] is around 0.4 for the present 
case p = 14/15) the layer is superlubric, i.e., there is no static friction (other than that, 
negligible, due to the weak 14/15 commensurability), and a time dependent modulation 
with amplitude Aq such that Vo +Aq < Vocut would produce no Shapiro steps. Thus we 
concentrate on corrugation amplitudes above the threshold, where the particle lattice is 
strongly pinned. 

Figure [H] is useful to understand the corrugation dependence of Shapiro steps. 
We probe three corrugation values: Vq = 0.5, Vq = 1 and Vq = 2 all with an 
amplitude modulation Aq = 70%Vo. The whole plateau structure moves to larger 
forces when the corrugation Vq is increased. As expected from common sense, the 
depinning force increases systematically with growing corrugation Vq, nearly ~ 20 times 
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Figure 11: p = 14/15: rescaled CM velocity for several corrugation amplitudes Vq, 
all with a modulation amplitude Aq = 70% Vq and frequency v = 0.02. The initial 
depinning force F increases rapidly with the corrugation amplitude Vq. Shapiro steps 
widen and shift to larger forces when the corrugation is increased. 


larger when Vq changes from 0.5 to 2. With reference to Fig. [5], both Vq = 0.5 and 
Vo = 1 have Vq — Aq < Vocrit, thus depinning occurs at an arbitrarily small force, 
also depending on frequency as discussed below in Sect. [5l In contrast, Vq = 2 has 
Vo — Aq = 0.6 > Vocrit, with a broad range of forces for which, independently of 
frequency, the system remains pinned, thus explaining the huge increase of static friction. 
The width of the Shapiro steps, especially the smaller ones, also increases with Vq: for 
example, the w/3 step, barely detectable for Vq = 0.5 and 1, becomes much broader 
for Vo = 2. In general, a larger potential Vq amplitude, corresponding to a comparably 
softer colloid lattice, spatially narrower solitons, and larger static friction, enhances the 
Shapiro-step structure, which emerges more clearly. 

Next, we explore the AC modulation amplitude dependence. Fig. [121 adopting 
Vo = 1 and varying Ao from 40% Vq, to Ao = 100% Vq (the purple triangles are indeed the 
same data as the purple triangles of Fig. [TT]). As in the lattice-matched case Fig. |3l when 
Ao increases, the Shapiro steps become wider and extend further down to smaller force 
F. Moreover, a larger AC modulation Ao depresses large-n longer-period subharmonic 
w/n in favor of shorter ones with smaller n (e.g w/2 or w). Conversely, as Fig. [12] 
shows, all individual plateaus shrink and become more numerous for smaller modulation 
amplitudes. Moreover, with Ao = 40% Vq we detect, in addition to w, w/2, w/3, and 
w/3, extra fractional plateaus such as Aw/b, 2w/3. In short, a small AC modulation 
amplitude Ao/Vq produces many weak subharmonic steps, while a large modulation 
promotes fewer stronger and wider ones. 

4.2. p = 3/4 

4 . 2 . 1 . Step structure. There is a clear contrast between the ’’nearly incommensurate” 
case, such as p = 14/15, or w = 1/14, and a ’’strongly commensurate” one, such as 
p = 3/4, corresponding to w = 1/3, Eq. (ITOll . In this strongly-commensurate case the 
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Figure 12: Velocity of the colloid CM (rescaled as in Fig. [2]), for p = 14/15, T = 0, 
oscillation frequency v = 0.02, corrugation amplitude Vq = 1, for different AC oscillation 
amplitudes: Aq = 40% Vq (dots), Aq = 70% Vq (triangles), Aq = 100% Vq (squares). 
Note that the smaller modulation Aq = 40% Vq exhibits a richer pattern of steps. The 
width of the main subharmonic steps tends to increase with increasing modulation 
amplitude Aq, at the expense of minor steps. 



Figure 13: p = 3/4, w = 1/3, N = 450 particles, corrugation Vq = 1, modulation 
Aq = 60% Vo, oscillation frequency n = 0.02. Two subharmonic Shapiro steps are 
visible. Inset (a): Magnihed view of the subharmonic plateaus. Inset (b): Magnihcation 
of the “trivial” plateau, Ucm/ (^^ oias) = 1 representing all colloids advancing as one every 
modulation period. 


Shapiro-step structure - Fig. [13]- shows fewer steps than in the nearly incommensurate 
case of Fig. |6l Basically only two subharmonic plateaus {w and w/2) plus the “trivial” 
step are detected. Moreover, the width of these plateaus is smaller than it was for 
p = 14/15. However, as was the case there, even for strong commensurability p = 3/4 
the plateau widths increase for increasing corrugation amplitude, and that even more 
rapidly than for p = 14/15. 
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Figure 14: A stroboscopic picture of the p = 3/4 advancement pattern, similar to 
Fig. [ini for the w/2 step. The antisoliton pattern is so weak that individual hexagonal 
domains are nearly invisible. The advancing colloids (bolder particles) are highlighted: 
150 particles advance in two sets at successive modulation periods: alternately 50 (left), 
then 100 (center), then again 50 (right), and so on. 


4-2.2. Subharmonic advancement of antisolitons. We examine the advancement 
mechanism for the subharmonic plateaus in Figure [TH The mechanism is again similar 
to that of Sec. 14.11 Fig. [101 Since p = 3/4 is further away from unity than 14/15, the 
antisoliton pattern is composed by much closer, partly overlapping and less distinct lines. 
The subdivision in domains is nearly invisible, with each particle always approximately 
as close to the edge of a domain as to its center. As a result, the oscillations around 
the average colloid trajectory are far less evident than in Fig. [HI Nevertheless, the 
advancement mechanics is the same as in the earlier case: e.g. the w/n subharmonic 
step involves Nw = 150 particles, advancing during n = 2 periods: we identify a clear 
pattern of advancing particles, displayed in Fig. [TH Here 150 colloids at the domain 
boundary advance in a first group of 50 particles in one modulation period, then a 
second group of 100 in the next period. In the w step, all 150 particles at the boundary 
cross the antisoliton line at every period: they could be identihed by superimposing the 
highlighted colloids in the hrst and second panels of Fig. [TH 

4 . 2 . 3 . Thermal effects. Finally we explore thermal effects on the subharmonic step 
structure. All tests in the p = 14/15 system show that at finite temperature the 
fragile Shapiro structure is strongly affected: the steps are washed away by the colloidal 
diffusive motion. No clear step can be detected for T > 0.001 -C T-^oom = 0.04, whether 
the standard Vq = 1 and Aq = 60% Vq or other parameters are adopted. 

By contrast, with the p = 3/4 mismatch ratio, where commensurability is stronger, 
the effects of thermal fluctuations are less dramatic, and at least the strongest steps 
survive. As shown in Fig. [151 fhe w plateau persists up to and above Tj-oom, the thermal 
effects leading to a slight reduction in width and to a slight edge rounding compared 
to T = 0. We explored even higher temperature, and found that even at T = dTroom 
the w step leaves some remnant of its existence in the force-velocity diagram. Fig. [T5h . 
However, even room-temperature fluctuations are enough to cancel out completely the 
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Figure 15: Thermal effects in the p = 3/4 mismatched conhguration, for Vq = 1, 
Aq = 60% Vq. Like in Fig. 01 the mean value and error bar of Vcm are computed by 
averaging over several simulations. The w step remains quite well identihed at Troom, 
while the w/2 step disappears altogether, (a): a zoom-in of the subharmonic w step at 
T = Troom compared with the much larger temperature T = 4Troom where the plateau 
structure is washed out by thermal noise. 
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Figure 16: Size scaling of the fluctuations of the CM velocity in room-temperature 
simulations of the p = 14/15 mismatched conhguration for Vq = 1, Aq = 60% Vq, 
executed for F = 0.08 (climb between steps) and F = 0.145 (step w) of Fig. [151 The 
various sizes are simulated with rectangular cells similar to the one depicted in Fig. [H 
The huctuations are obtained as the standard deviation evaluated over 4 independent 
simulations for each size. Dashed line: a guide to the eye. 
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w/2 plateau. 

Thermal effects on Shapiro steps also depends to some extent on size. As long 
as the colloid density remains uniform, thermal huctuations average out for increasing 
sample size and averaging time: for a given temperature, a larger sample size and longer 
simulation show hatter, sharper plateaus. Figure [16] shows how the thermal huctuations 
on the colloid CM speed depend on the number N of colloids in the simulation supercell. 
For large size, as expected by statistics, thermal huctuations scale as while when 















Subharmonic Shapiro steps of sliding colloidal monolayers in optical lattices 


22 



n ^ . ^ 




. • • •' 

c 


t 

fc. 

(U 

M'' ^ . 

Ui 

O - . 


tL- 0.1 


Oh 

c 



\ 

\ 

• 

*o 

II 

c 

■q- 

♦ ■■■ v” 

Q 

□ p=14/15, w step 


♦ p=14/15, w/2 step 

0.01 

A p=14/15, w/3 step 


•a V 


1 , , .1 , , , , 


0.001 0.01 
Frequency v 


Figure 17: The static depinning force as a function of the modulation frequency in 
conditions where the sliding barrier vanishes during the modulation period. Circles: the 
lattice-matched geometry, with Vq = 1, Aq = 100% Vq. Other symbols: p = 14/15, with 

K, = l, Ao = 70%K). 


the cell lateral size is comparable to or smaller than the correlation length associated 
to the colloid-colloid interaction, dynamical correlations enhance the fluctuations. Note 
that inside a plateau (triangles in Fig. fTOjl . the fluctuations are much smaller and the 
correlation length (marked by the knee in the size scaling) signihcantly longer, than in 
a climb between steps. 


5. The depinning force 


This section deals with the DC-Force - corrugation-amplitude diagram. Fig. [5l We 
can foresee that, depending on Vq and Aq, the static-friction force Fhepm, namely the 
minimum force required to maintain a net nonzero sliding speed, can fall in one of two 
alternative regimes. The simplest condition occurs when Vq — Aq exceeds the Aubry 
Vocrit^ then the colloid remains statically pinned throughout the modulation period, with 
the result that sliding requires a hnite driving force greater or equal than the depinning 
force appropriate to a corrugation Vq — Aq in the absence of modulation. 

Conversely, when Vq — Aq < Vocrit, for a fraction /^ee of the oscillation period 
the pinning barrier vanishes, and the colloid slides freely. A slow enough AC modulation 
gives the particles enough time to take advantage of the temporary lack of pinning and 
generate the hrst Shapiro step for arbitrarily small driving force F. In such condition, 
the depinning force is a function of the modulation frequency. The depinning force Fdepin 
can be estimated as follows: During the unpinned time /free^^”^ the colloids advance at 
a speed ~ pF, where p ~ rj~^ is the colloid mobility; by equating the displacement 
fireei^~^Fp to a minimum advancement compatible with the overall displacement in the 
lowest Shapiro step, say ^uias, we obtain 


depin — 


^®las 

^/^/free 


(14) 
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Thus, under modulation conditions such that Vq — Aq < Vocrit, the static friction force 
is predicted to vanish linearly with the modulation frequency. Simulations conhrm this 
power law scaling, indicated by the dashed line in Fig. [171 For decreasing frequency, 
simulations also show that the n > 1 subharmonic steps tend to weaken and eventually 
disappear, see e.g. the diamonds (w/2) and triangles (tc/3) in Fig. UTl As a result, 
eventually at very low frequency depinning occurs directly into the w step. 

The fully-matched case p = I discussed in Sec. |3] has no Aubry transition, so 
hocrit = 0. Even so, at 100% modulation rate (Aq = Vq), in the brief instants when 
sin(ci;f) ~ —1, again free sliding is allowed. At hnite driving force F a hnite time interval 
/free^^”^ is Compatible with barrier overcoming and consequent free sliding. This interval 
is obtained by evaluating the region of F > F^i using the expression ([9]) (under the 
condition Aq = Vq) for fhe barrier crossing force, obtaining 

11 . f9Fai^ \ 

/free = Tt + “ ^rCSin - 1 . 15 

2 TT V SttFo / 

By expanding the arcsin function with its argument near —1 (small driving force) we 
obtain /free — [OFoias/(dvr^Vo)]^/^. As above, we then equate the displacement ffreeJ^~^Fp 
to the minimum advancement required for the colloids to advance by one corrugation 
lattice spacing, say |aias, obtaining 

iAepin ^ TT • (16) 

Again, simulations conhrm this result, see Fig. [T71 with depinning always occurring into 
the trivial harmonic step. In principle a similar “critical” power-law regime is to 
be expected even in the mismatched case p < 1, under the condition Vq — Aq = Vocrit 
of an instantaneous touching of the free-sliding superlubric Aubry phase during each 
modulation cycle. 

6. Discussion and Conclusion 

The present work addresses, by means of MD simulations, the question whether Shapiro 
step structures can occur in the forced sliding of a 2D colloidal monolayer immersed in an 
optical lattice spatially periodic corrugation potential whose amplitude has a static DC 
component plus an AC modulation that oscillates sinusoidally in time. We conhrm that 
Shapiro steps should appear, in the form of intervals of static force within which the CM 
colloid velocity is quantized, at least at zero temperature, to a value entirely determined 
by the modulation frequency and by geometric parameters. The only condition required 
for this to occur at T = 0 is that in the absence of the AC modulation the colloids should 
be pinned by static friction, and not free sliding, as can happen when the 2D colloid 
lattice and the optical lattice are mismatched and the optical lattice is sufficiently strong. 
In commensurate, perfectly lattice-matched conditions p = 1 only ’’trivial” integer steps 
are expected, corresponding to the synchronized periodic advancement of all particles at 
each AC modulation period. In lattice-mismatched conditions p % 1 we find additional 
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subharmonic Shapiro plateaus beside the integer ones. The crucial quantity governing 
the subharmonic plateaus is w = |1 — which determines the quantized velocity step 
value in the totally general form 

Vcm = ai^U! — , (17) 

n 

where n is the AC frequency, oias is the optical lattice spacing, and m and n are small 
integers. Interestingly, we hnd that every quantized Shapiro step has a specihc real- 
space signature corresponding to which particles advance synchronously with the AC 
modulation - a signature which is purely geometric. In a step the total number of 
advancing particles per modulation period equals exactly a fraction Nwm/n of the 
total particle number N. However, a step corresponding to integer n lasts n periods. 
During each period only a sub-fraction of these particles advance by one lattice 
step, and the total Nwm is achieved only at the end of a cycle of n periods. 

We analyzed how the corrugation amplitude Vq and modulation Aq influence the 
ensuing step structure. An increase in Vq leads to a comparably softer colloidal layer, 
thus increased localization and static friction, resulting in a shift to higher forces F of 
the Shapiro plateaus and generally an increase in their width. On the other hand an 
increase in Aq to the point of making Vq — Aq < Vocrit leads the oscillation to explore 
more and more of the hard-colloidal-layer region for part of the oscillation period. In 
this region, free sliding is possible, as well know in the ID FK model [3H]. This results 
in more extended steps, in particular down to smaller F. In such cases, the smallest 
(depinning) force is a function of the modulation frequency, and we verihed that it 
changes either as u as depending on whether a free sliding state is crossed for a 
hnite fraction of the modulation period or just touched instantly. Additionally, a very 
large Aq — Vq tends to simplify the rich ladder of subharmonic step structures observed 
at intermediate Aq. 

We also checked how important the specihc shape of the corrugation-amplitude time 
modulation really is for the possibility to observe Shapiro steps. We ran simulations 
both for the lattice-matched case and for the p = 3/4 mismatch ratio, replacing sin(a;t) 
in Eq. ((Tj) with a square wave sign(sin(a;t)), where sign(a;) equals x/\x\ for nonzero 
X. Compared to the sinusoidal case with identical modulation amplitude Aq and 
other conditions, the square-wave simulations show quite similar velocity-force prohles, 
with marginally wider and more robust Shapiro steps. This small difference can be 
understood as a consequence of the larger RMS modulation of the square wave than the 
sinusoidal wave for the same peak-to-peak amplitude. 

Thermal huctuations and density inhomogeneities present in real colloid sliding 
experiments are likely to affect the possibility to observe Shapiro structures. While 
the perfectly-matched system is barely affected by the thermal Brownian motion, the 
subharmonic Shapiro-step structures of mismatched systems are far more fragile. Near 
matching between the colloid and the optical lattice, the loose soliton network is strongly 
affected by thermal fluctuations, which wash easily away most of the delicate Shapiro 
plateaus. By contrast, for a more robust (but commensurate) lattice mismatch such as 
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p = 3/4, we retrieve subharmonic steps up to Troom and even above. We test also the 
size dependence in the latter case: as long as the density is uniform, a greater number 
of particles averages out thermal fluctuation, ending up in a flatter and less rounded 
step. 

The model explored in the present work neglects hydrodynamic interactions [53] : 
while such velocity-dependent many-body forces may affect the quantitative detail of 
the Shapiro steps, we do not expect they would change the qualitative picture much. In 
particular, nonuniform solitonic motion, as occurs in mismatched conhgurations, could 
be enhanced by hydrodynamic interactions. We are currently investigating this issue in 
quantitative detail. Also, the detail of the two-body interaction, Eq. ([5]), is not likely 
to play any qualitative role in the Shapiro-steps physics. In particular, we expect a 
similar behavior if the colloid-colloid interaction, rather than screened electrostatic, was 
the power law describing the repulsion between magnetic colloidal particles in a layer, 
controlled by a perpendicular magnetic field [Ml ESI [56] • 

In conclusion, this theoretical study predicts that Shapiro-like plateau structures of 
colloid CM velocities should be easily detected experimentally in a 2D rich geometry- at 
least in lattice-matched condition, with potential amplitude and oscillation modulation 
values varying over broad ranges. In mismatched conditions, on the other hand, the 
experimental detection of subharmonic plateaus seems limited by the ability to provide 
large regions where the 2D colloid density is constant in space with sufficient accuracy. 
While experimentally this may not be possible for all subharmonic quantizations, still 
for a lattice-spacing ratio sufficiently distant from unity, a step structure should be 
detectable under experimental conditions. Subharmonic steps were indeed detected 
(although not identihed as such) in the simpler ID experimental setup of Ref. [H] . 
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